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Abstract 

X-ray scattering experiments on bulk liquid 3 He fl]-f|] have indicated the possibility of the ex- 
istence of a sharp collective mode at large momentum transfers. We address this issue within a 
manifestly microscopic theory of excitations in a Fermi fluid that can be understood as proper 
generalization of the time-honored theory of Jackson, Feenberg, and Campbell of excitations 
in 4 He. We show that both neutron and X-ray data can be well explained within a theory where 
the high momentum excitations lie in fact inside the particle-hole continuum. "Pair fluctuations" 
contribute a sharpening of the mode compared to the random phase approximation (RPA) . When 
the theoretical results are convoluted with the experimental resolution, the agreement between 
theory and X-ray data is quite good. 

PACS numbers: 67.30.-n, 67.30.em 
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I. INTRODUCTION 



The helium fluids 3 He and 4 He are the prime examples of strongly correlated quantum 
many-body systems. Governed by a simple Hamiltonian, yet very dense, they have been 
studied for decades and still offer surprises leading to new insights. It is fair to say that 
understanding the helium fluids lies at the core of understanding other strongly correlated 
systems. 

Liquid 3 He is for many reasons the more challenging substance for both, theoretical and 
experimental investigations. Key information on the dynamics of the system is provided 
by the dynamic structure function S(q;uj). Experimentally, the dynamic structure function 
of 3 He is mostly determined by neutron scattering. The results are well documented in a 
book f?J, the theoretical and experimental understanding a decade ago has been summarized 
in Ref. |[ Recent inelastic X-ray scattering experiments have led to a controversy on the 
evolution of the zero sound mode at intermediate wave- vectors [fj]-|3|. It is one of the purposes 
of this paper to examine this issue on the basis of microscopic many body theory. 

The ground state of liquid 4 He is well understood from both diagrammatic many-body 
techniques and simulation methods, see Ref. |9| for a collection of pedagogical articles. Build- 
ing upon pioneering work by Jackson, Feenberg, and Campbell recent developments 
|T0|^T2[ have brought manifestly microscopic theories of 4 He to a level where quantitative 
predictions for the excitation spectrum are possible up to momentum transfers well beyond 
the roton minimum. 

We have recently generalized the theory of Refs. [|-[] to Fermions [IB]. From the the- 



oretical point of view, the major technical problem is to deal with the multitude of ex- 
change effects in an efficient manner; we have formulated a tractable strategy for including 
pair-excitation effects in Fermi fluids. This paper is concerned with a theoretical study of 
excitations and damping in normal liquid 3 He at atomic wave lengths based on that theory. 

II. DYNAMIC MANY-BODY THEORY 
A. Equations of Motion 

This section gives a very brief compilation of the techniques developed in Ref. [13]. We de- 
liberately refrain from any explanation and refer to the original work and earlier pedagogical 
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material || for details. 

Microscopic many-body theory starts from a phenomenological Hamiltonian for N inter- 
acting fermions, 

^ = -E^V J 2 + E-(|r 4 -r J |) • (2.1) 

» *<?' 

For strong interactions, correlated basis functions (CBF) theory || has proved to be an effi- 
cient and accurate method for obtaining ground state properties. It starts with a variational 
wave function of the form 

F 1$ ) 

m ) = 1 o/ (22) 

where $ (1, ■ ■ ■ ,i, ■ ■ ■ , N) is a model state, normally a Slater-determinant, and "z" is short 
for both spatial and v discrete (spin and/or isospin) degrees of freedom. The correlation 
operator F(l, . . . , N) is suitably chosen to describe the important features of the interacting 
system. Most practical and successful is the Jastrow-Feenberg || form 



F(l,...,N) = exp<{i 



l<i<j<N l<i<j<k<N 



(2.3) 



The u n ( r ij • • ■ j r n) are made unique by requiring them to vanish for [r, — r 3 -| — t-oo ("cluster 
property"). 

From the wave function ( |2.2|) , ( |2.3|) , the energy expectation value 



H 0>0 = <* |#|*o) (2.4) 

can be calculated either by simulation or by integral equation methods. The hierarchy 
of Fermi-Hypernetted-Chain (FHNC) approximations is compatible with the optimization 
problem, i.e. with determining the optimal correlation functions u n (ri, . . . , r n ) through 
functionally minimizing the energy 

6H °'° -0. (2.5) 



5u n (rx, . . . ,r n ) 

The dynamics of the system is treated with the logical generalization of the wave function 
Q2.2|) , ( |2.3|) by writing the response of the system to a weak and time dependent external 
field H ext (t) = Ei 5/i ex t(ri; t) in the form 



|* (t)> = FeW*)|$ > , ^ = <* (t)|* (t)>. (2.6) 



where F is taken from the ground state calculation, and SU(t) is a sum of n-particle-n-hole 
excitation operators: 



6U (*) = 5u P h{t) a l a h + | 6u PP 'hh'(t) a\a\,a h ,a h + . . . ; 

ph pp'hh' 



(2.7) 



As a convention, we will label "particle" (unoccupied) and "hole" (occupied) states with 
p, p' and h, h', respectively. Discrete degrees of freedom are suppressed. 

Equations of motion for the 5u^\t) are derived from the linearized action principle [0, |TB 

d 



5 JdtC{t) = 5 JdtU>{t) 



H + H ext (t) -m 



Of 



9(t) 







(2.8) 



The analytic and diagrammatic manipulations to bring the resulting equations of motion 



into a numerically tractable form are lengthy and delicate |13[ , but the result is surprisingly 
simple and can be cast into the familiar form of the time-dependent Hartree-Fock (TDHF) 
theory |16|] with energy dependent effective interactions. The induced density fluctuations is 
expanded in terms of matrix elements of the density operator in the non-interacting system, 
(h\5p(v)\p). 



Po,ph( r 



5p(r; w) = £ £ [pJU( r ) S<^\u) + Po* ph (r) 5c ph \u) 

ph 



(2.9) 



The amplitudes 5c^\oj) are related by the THDF equations to the matrix elements hi 
(i\6h ext \j) of the external field. 



(hu - e ph )d p , p >d h , h > + V p { fy h ,(uj 



Vpp'hh'o( U ) 



6a 



v 0, pp'hh' 



-(hu + e ph )6 PtP ,6 h , h/ + V ( p $ tph (uj) j I 8c K hp 



(+) 

ph 

,(-) 



hph 

hhp 
(2.10) 



where e p h = t(p) — t(h) is the particle-hole excitation energy of the non-interacting system 
and t(q) = h 2 q 2 /2m is the kinetic energy. 

Thus, the machinery of microscopic many-body theory has led to a definition of the 
effective interactions V ph ^ p , h ,{uj) and V ppl ^ hh ,(uj) in terms of the correlation operator F and 
the underlying microscopic Hamiltonian. 



B. Local interactions 



Much of our qualitative understanding of the dynamics of interacting quantum fluids is 
based on the "random phase approximation" (RPA). In that approximation, the density- 
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density response function is of the form 

X { q -u)= ~ Mq ' U) (2.H) 

i - U?) Xo(g; w) 

Here, XoC?;^) is the Lindhard function, and %- h (q) an appropriately defined static "particle- 
hole interaction" or "pseudo-potential [fL7|'\ ^ * s related to the dynamic structure function 



S{q;u) = --$Sm\x{q',u>)]e(u) (2.12) 

7T 



which satisfies, amongst others, the sum rules 



m (q) = S(q) = / dhu S(q;u), (2-13) 
J o 

m i(q) — — / dhu hu S(q;u) , (2-14) 

2m jo 

where S(q) is the static structure factor. A version of the RPA with an effective inter- 
action can be derived by restricting the excitation operator ( |2.7[ ) to one-particle-one hole 
amplitudes, and omitting (or approximating) exchange terms. 

The RPA (]2 . 1 1|) for x(q] w ) displays the essential features of the dynamic structure func- 
tion S(q; u) qualitatively correctly: S(q; u) can be characterized as being a superposition 
of a collective mode similar to the phonon-maxon-roton in 4 He, plus an incoherent particle- 



hole band which strongly dampens this mode |L8[ where this is kinematically allowed. A 
schematic drawing is shown in Fig. [l|. 

The picture misses important physics: In 3 He the RPA, when defined through the form 
( P- 11| ) and the sum rules ( 2.13 )— ( |2.14|) , predicts a zero-sound mode that is significantly 



too high. This is consistent with the same deficiency of the Feynman spectrum e(q) = 
h 2 q 2 /2mS(q) in 4 He. Drawing on the analogy to 4 He ||18|| , this deficiency can to a large 



extent be cured by introducing pair fluctuations u^, hh ,(t) in the excitation operator. 

In an approach that maintains the locality of the effective interactions entering the TDHF 
equations, we have in Ref . [13] derived practical working formulas for the effective interactions 
in terms of a three-body vertex and a pair propagator, these are immediately recognized as 
generalizations of the bosonic version: 



N 

V p %» = 6p+P ' N h - h ' V B (q;co). (2.15) 
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FIG. 1: (color online) A schematic picture of the dynamic structure function S(q;uj) in 3 He. The 
gray shaded area is the particle-hole continuum of non-interacting particles, the circles show the 
maximum intensity of neutron-scattering data (Ref. |l~9| ) and the dots are X-ray scattering data 



1. The red line shows the results of our calculation. See also Sec. Ill for further discussions. 



The components of the (energy-dependent) interactions V AB (q;uj) are 



VMM 



^p-h(g) 

V p - h (q) 



[*}) 2 W A (q;u) + [a;) 2 W*(q; -co) 
W A (q;u) + W:(q; -u) 



(2.16) 
(2.17) 



with = [S F (q) ± S(q))/2S(q). Here, the V^,_h(?) is the static part of the particle-hole 
interaction that is related to the static structure function through the RPA relationship 
( P- 11| ) and the two sum rules ( |2.13 ) and (|2.14 ). The energy dependent part of the interaction, 
W A (q; u) describes the splitting and re-combination of phonons, it consists of a three-phonon 
vertex K q , q > q " and a two-phonon propagator £ ,_1 (g / , q"; uj). The energy dependent correction 
to the effective interaction is 



(2.18) 



with a three-body vertex 



Kq,q'q" 



h 2 



S(q')S(q" 



2m S F (q)S F (q')S F (q") 



q- q'l dd (</) + q- q"X dd (g") - q 2 u 3 (q, q' , q")} .(2.19) 
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Here, S F (q) is the static structure function of non-interacting fermions. X d d(q) = ^-/S F (q) — 
1/S(q) can be identified with the set of "non- nodal" diagrams of the Fermi-hypernetted- 



chain theory, and Us(q,q',q") is the three-body ground state correlation [|20 
Finally we need the pair propagator: 

-— - re(<7i;u/) n(q 2 ;uj-uj') (2.20) 
2m 

— oo 

K(q;oJ) = ~ (2.21) 

1 + huT dd (q)n (q; u) 

with the "direct-direct" correlation function rdd(?) — (S(q) — S F (q))/S F (q), and the partial 
Lindhard function: 

K (q; uj) = 1 V - npHh . . (2.22) 

ovy ' J N ^ hu - e ph + ir] K J 

Evidently, we need for the execution of the theory only the static structure function S(q). 



In this work, we have used input from the FHNC-EL calculations of Ref. ^ that were based 
on the Aziz-II potential ||21|| . 

C. Bose limit 

For the sake of discussions and to make the connections with previous work we briefly 
discuss the Bose limit. In that case, S F (q) = 1, and we can identify rdd(qO —±S(q) — l. The 
constituents of the pair propagator become 

1 



C se (^) 



hco — t(q) + if] 



K Bose (q;u) = -i— ■ (2.23) 

The pair propagator is then simply 

^iose(9i, <l2] w) = 1 — . . (2.24) 

S(qi)S(q 2 ) hu - e(qx) - e(q 2 ) + it) 

The form (p.24 ) shows more clearly than Eq. ( p.20[) the physical meaning of the pair 



propagator. It also shows the consequence of restricting the excitation operator to pair 
fluctuations: The pair propagator contains in the energy denominator the Feynman spectrum 
which is known to be almost a factor of two higher than the physical phonon-roton spectrum. 



For bosons, the problem has been cured by introducing multi-particle fluctuations Jill 12 



these investigations show that the energy denominator gets indeed renormalized and closer to 
the physical spectrum. We have in the past estimated the effect of higher-order fluctuations 
by scaling the Feynman spectrum in the energy denominator such that the energy of the 
scaled spectrum agrees roughly with the experimental roton minimum, i.e. we have used 

111 

E\ (qi,Q2;u) = —— . = -E~ (q 1 ,q 2 ;u)/\) . (2.25) 

S{qi)S{q 2 ) tiuj - Ae(gi) - Xe(q 2 ) +irj X 

We will use this procedure below to estimate the effect of multi-particle fluctuations in 3 He. 



III. DYNAMIC STRUCTURE OF 3 HE 

We have carried out a comprehensive array of calculations of S(q; u) for 3 He in the exper- 
imentally accessible density regime. Fig. |1] shows, at saturation density, the collective mode 
and the particle-hole continuum corresponding to single-particle energies of non-interacting 



particles. Our theoretical results are compared with X-ray |1]] and neutron scattering 
data. The "zero sound" mode is, when inside the particle-hole continuum, identified as the 
maximum value of S(q; u>). 

Figs. | give a first glance of our results. The figures show S(q; to) for the densities 
p = 0.0166 and 0.020 A -3 , this corresponds to SVP and at a pressure of 10 bar, respectively 



22] , [2311 . The figures show the typical scenario of a collective mode at long wave lengths, 
and a particle-hole band. Note, however, that our particle-hole band does not have sharp 
boundaries, this effect is due to pair fluctuations which have the consequence that the 
interactions V AB (q;uj) can be complex. 

Recent X-ray scattering results [|1] have led to a discussion on the location of the particle- 
hole band 0, |3[]. The authors of Ref. p] state that The obtained results show no evidence of 
such a decay: the zero-sound mode remains well defined in the whole explored wave number 
range. According to that analysis, the picture shown in Fig. [i] would not be correct, 
rather the particle-hole band should be significantly lower in energy. In particular the 
collective mode would lie outside the particle- hole band (See Fig. 3 of Ref. [I]). This is 
somewhat surprising because such a scenario would require a single particle spectrum that is 
characterized by an average effective mass m* w 3m up to more than three times the Fermi 
momentum. This is in contrast to theoretical work that predicts that the effective mass 
should have a peak around the Fermi momentum caused by spin-fluctuations, a secondary 
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FIG. 2: (Color online) -9mx(?;w) at a density of p = 0.0166A~ 3 (left pane) and p = 0.02A" 3 
(right pane). The contour lines are at 0.02, 0.05, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0 tp 1 
where tp = fi 2 kp/2m is the Fermi energy of non-interacting 3 He atoms. 

maximum at about twice the Fermi momentum due to density fluctuations, but then falls 



off rapidly towards the bare mass [24-27 



Further details are shown in Figs. |3|. To facilitate the comparison with experiments, we 
have convoluted our theoretical spectra with the experimental resolution. We have scaled 
our results by 1/S(q) such that the integrated strength is 1 for all momentum transfers. We 
have also scaled the experimental spectra, the scaling factor agrees, apart from very long 
wave lengths, with our 1/S(q). 

Evidently, our results agree, after introducing the broadening, rather well with the X-ray 
scattering data. It appears therefore that the X-ray data can be explained equally well by a 
scenario where the collective strength is inside the particle-hole band. We stress, however, 
that such a good agreement can not be obtained in an RPA model that assumes a form 
( [2. 11| ) for the density-density response function and determines the effective interaction by 
the sum rules ( |2.13|) and ( |2.14|) . Typically, the strength is distributed over a wider regime, 
thus suggesting a broader maximum than observed. 

Fig. |3|b shows the same results for the density of p = 0.02fjA -3 . 
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FIG. 3: The left pane shows the normalized dynamic structure function S (q; lo) / S (q) of 3 He at a 
density of p = 0.0166A" 3 , for a sequence of wave numbers. Our theoretical results (solid line) are 
compared with the X-ray scattering data 0] (squares) at saturated vapor pressure. Also shown 
are our theoretical results convoluted with the experimental resolution (dashed). The right pane 
shows the theory at an density of p = 0.02A -3 . All spectra are normalized to unity. The secondary 
maximum observed in Fig. || at shorter wave lengths is the due to the possible decay of excitations 
into two "maxon"s and will likely be broadened by multi-pair fluctuations. 

A. Multi-pair coupling 



The question arises naturally whether there could be an equivalent to the "Pitaevskii 
plateau" |?8 in 3 He similar to one in 4 He. In 4 He, this plateau is caused by the kinematic 
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possibility that a "phonon" can decay into two "roton" excitations. The same process can, 
of course, also happen for fermions. However, one expects that the resulting strength is 
smeared out because the "roton" is per se not a sharp excitation but can be observed, at its 
best, by a high density of states in the area where the roton would be in a Bose system. 

To determine if such an effect could be visible we recall our discussion of the Bose limit 
in Sec. |11 (J| , in particular the the consequences of restricting the excitation operator ( |2.7|) to 



pair fluctuations. The consequence of this version of the theory is that, in a Bose fluid, the 
"plateau" appears at an energy hu p \ = 2e(qA) where is the wave number of the "Feynman 
roton", i.e. the minimum of the Feynman dispersion spectrum e(q) [fLOR . We should therefore 
also expect such a feature for fermions only at about the energy of the maximum density 
of states in the RPA. In both the RPA response function and in the propagator n(q; oS) this 
energy is of the order of 1.4 meV. The position of the "roton" is lowered by pair-fluctuations 
to less than 1 meV, cf. Fig. |[ 

In order to estimate the consequence of including multi-particle fluctuations, we have 



in previous work [^9| scaled the Feynman spectrum in the pair propagator such that the 
spectrum in the energy denominator roughly agrees with the result of the CBF calculation. 
Without such a scaling, one would, in the present calculation, expect a plateau-like structure 
at an energy of hu p \ ~ 2.8 meV and not at twice the energy of the actual "roton". We adopt 
here the same philosophy as in the work on bosons: We scale the two-phonon propagator 
as formulated in Eq. ( |2.25| ) . This modification has little effect for long wave lengths, but it 
changes the features of S(q; u>) for short wave lengths visibly as seen in Fig. f|. The most 
notable effect is the expected broadening of S(q; uS) at high momentum transfers. 

A closeup of the situation is given in Fig. |5] where we show S(q; u) at a relatively large 
momentum transfer q = 4k f- Obviously the more realistic description of two pair excitations 
provided by the scaling procedure ( |2.25|) has the effect of smoothing the spectrum and, in 
particular, adding strength to the low-energy part of the spectrum which improves the 
agreement between experiment and theory further. 

IV. SUMMARY 

We have in this paper presented the application of a new theory of the dynamic structure 
function of Fermi systems to 3 He. The essential results have been discussed in the text. Let 
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FIG. 4: (Color online) Same as Fig. || where the pair propagator is scaled by a factor A = 0.65 as 
described in Eq. ( |2.25| ) 




FIG. 5: (color online) The figure shows the dynamic structure function of 3 He at saturated vapor 
pressure at q = 4kF- The solid blue line is the result of the pair excitations theory and dashed 
red the RPA result. The magenta line shows the result for the pair excitations where the pair 
propagator has been modified by a factor A = 0.65 as described in Eq. ( |2.25 ). The corresponding 
dashed lines are the theoretical result folded with the experimental resolution. Inelastic X-ray 
diffraction data obtained by Albergamo et al. Q are also shown (boxes). 



12 



us reiterate here first the basic theoretical objectives of our work: The key to a quantita- 
tively correct description of the dynamics of 3 He - and, therefore, other strongly correlated 
quantum many-body systems - is the fact that, for excitations as atomic wave lengths, the 
short-ranged structure of the wave function is time dependent. The intermediate states can 
be described only in terms of the quantum numbers of at least two particles. With that 
strategy we have achieved a rather satisfactory agreement between experiment and theory. 

Our results have been obtained without any modification of the particle-hole spectrum. Of 
course, we do not claim that the precise location of the single-particle spectrum is completely 
irrelevant for the energetics of the zero sound mode; but we have demonstrated that the 
dominant mechanism in Bose and Fermi fluids is the same, namely pair-fluctuations. 

In independent work [^, |27|, [3(J, we have used the ideas of CBF theory as well as the 
Aldrich-Pines pseudopotential theory to calculate the single-particle propagator in 3 He. In 
both three and two dimensions, we found the above-mentioned strong variation of the ef- 
fective mass with momentum. Near the Fermi surface, we found good agreement between 
the theoretical effective mass and the one obtained experimentally from specific heat mea- 
surements. It is, however, not clear at all if the simple concept of a single particle spectrum 
with a possibly momentum-dependent effective mass is appropriate away from the Fermi 
surface. Even simple concepts like the G(0)W approximation leads to an energy depen- 
dent self-energy. In order to maintain the sum rules ( |2.13| )-(2.14), any modification of the 



particle-hole spectrum must go along with an inclusion of exchange effects. At the level of 



single-pair fluctuations such a calculation is quite feasible [pl| , p2| but its implementation at 
the level of the theory described above still needs to be done. 

In terms of the interpretation of experiments we hope that we have contributed to the 
clarification of X-ray scattering results. In view of the large experimental width one is 
inclined to consider the theoretical description as satisfactory. At longer wave lengths there 



are still some open questions with respect to the onset of Landau damping |19| which we 



hope to clarify in the future by the inclusion of the above-mentioned self-energy corrections. 
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